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We investigate static properties of laser-driven, ultracold Rydberg atoms confined to one- and 
two-dimensional uniform lattices in the limit of vanishing laser coupling. The spectral structure of 
square lattices is compared to those of linear chains and similarities as well as differences are pointed 
out. Furthermore, we employ a method based on elements of graph theory to numerically determine 
the laser detuning-dependent ground states of various lattice geometries. Ground states for chains 
as well as square and rectangular lattices are provided and discussed. 
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I. INTRODUCTION 

Recent years have seen an increasing interest in the 
coherent control of ultracold Rydberg states, i.e., states 
with large principle quantum number n with energies just 
below the ionization limit. The vast displacement of the 
valence electron leads to a large polarizability and there- 
fore makes Rydberg atoms very susceptible to external 
fields Furthermore, these polarizabilitics can lead 
to very strong Rydberg-Rydberg interactions of dipole- 
dipole or van-der-Waals character. These interactions 
can be tuned by applying external fields as well as by 
the choice of the atomic state. The so-called blockade 
mechanism resulting from the strong mutual interaction 
is an essential building block in many works Q, as it is 
in the present one. In a nutshell, the Rydberg blockade 
mechanism means that the strong Rydberg-Rydberg in- 
teraction induces a level shift of many-body states with 
more than one Rydberg atom present such that the ex- 
citation of an atom in the vicinity of a Rydberg atom is 
shifted off- resonant and is therefore prevented (blocked). 
The dipole-blockade mechanism has been predicted the- 
oretically in works addressing quantum information pro- 
cessing for the implementations of quantum logic gates 
@, Hi- Few years later, various groups succeeded to 
demonstrate the blockade mechanism also experimentally 
in ensembles of rubidium atoms 043 and cesium atoms 
Q. More recently, two experiments demonstrated the 
Rydberg blockade effect for two individual atoms with a 
distance of several micrometers 0, [loj and subsequently 
the realization of a CNOTgate and the generation of en- 
tanglement [ll|, EH • In [13, EH the associated enhanced 
collective Rabi frequency has been experimentally veri- 
fied. 
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The ability to precisely control the trap ping potential 
of Rydberg atoms by means of magnetic [14j and opti- 
cal [15[ traps in combination with the long-range mutual 
interaction put a lot of attention of theoretical works 
to Rydberg atoms in lattices. For instance, in fl6l - fl8| 
a one-dimensional lattice with periodic boundary condi- 
tions, i.e., a ring lattice, is considered. In the regime 
of a weak laser coupling with respect to the interaction 
strength of neighboring Rydberg atoms, the time evolu- 
tion of Rydberg densities and correlations has been an- 
alyzed, whereas in the regime of a strong laser coupling 
an analytical approach lead to collective fermionic excita- 
tions. In a recent work we have presented spectral prop- 
erties of finite laser-driven lattices with uniform as well as 
variable spacings [l9| . This knowledge can be exploited 
to determine the strength of the Rydberg-Rydberg inter- 
action by means of the excitation dynamics poj . 

Several works investigate the ground states of Rydberg 
lattices since they are premier candidates for an exper- 
imental realization. In [2l[ a method to selectively ex- 
cite crystalline structures of Rydberg atoms by employ- 
ing a chirped laser pulse is proposed. The ground state 
of a one-dimensional lattice has been analyzed in (22j by 
treating the influence of a finite laser coupling perturba- 
tively. The resulting phase diagram of the system shows 
a quantum melting of the crystalline phases. There are 
also first theoretical works on ground states of two di- 
mensional lattices, which indicate phase transitions from 
an ordered to a disordered phase [231 ] . 

The present work expands our previous investigations 
of laser-driven lattices of ultracold Rydberg atoms to two 
dimensions and focuses on their ground states in the limit 
of vanishing laser coupling. Because of the enormous 
state space, sophisticated mathematical tools are nec- 
essary to handle the determination of the ground state 
even for the simple case of vanishing laser coupling. The 
basic model, the associated Hamiltonian and its symme- 
tries are provided in Scc.|TTJ Subsequently, in Sec. Mil we 
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compare the spectra of chains and square lattices in the 
limit of vanishing laser coupling. Section HVl which is ad- 
dressing the laser detuning-dcpcndent ground states, is 
divided as follows: After defining the ground state prob- 
lem, we show in Sec. IIV Al that it is equivalent to the so- 
called maximum cut problem. Section llVBI is dedicated 
to the numerical results of the maximum cut problem, 
wherein we provide the ground states for some chains, 
square lattices and an example of rectangular lattices. 
We conclude with a brief summary (Sec. [V]) and an ap- 
pendix (Sec.[X|), which contains a more detailed calcula- 
tion concerning the Ising-type form of the Hamiltonian 
used in Sec. [TVl 



II. THE MODEL 

Throughout this work we consider the same model as 
in [2(|: An ultracold gas is trapped in a finite, regu- 
lar one- or two-dimensional optical lattice with in total N 
sites, each containing Nq atoms. An additional external 
laser field couples the ground state to a Rydberg level by 
means of a two-photon transition with a Rabi frequency 
of Oo. If the excitation lasers are far detuned from the 
intermediate level, it can be adiabatically eliminated as 
demonstrated in |24| . We also allow for a small detuning 
of the two-photon transition with respect to the Rydberg 
level. Because of the adiabatic elimination of the inter- 
mediate level, the single-atom model so far reduces to 
a two-level system. When assuming in addition that the 
blockade radius is larger than the spatial extension of the 
individual sites, each site k can again be described as a 
two- level system with a collective ground state | g ) k and 
an excited state | e ) fc which is characterized by a symmet- 
rically shared Rydberg excitation. The coupling between 
the levels of this so-called superatom is thus given by the 
enhanced collective Rabi frequency 51 = ^/NqUq. 



A. Hamiltonian 

Taking into account the Rydberg interaction between 
different sites, the final model Hamiltonian of the system 
reads in the basis of the superatom states {\g) k , | e ) ^ } 
and after applying the rotating wave approximation 
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Here <Jx , ai are the Pauli matrices and rie 
1] is the excitation number operator, where the label k 
denotes the site index. Equation ([l) can also be inter- 
preted as the Hamiltonian of an interacting spin chain 
subject to two external orthogonal magnetic fields. 



The first term H^ct = y SfcLi a 'z'' contains the laser 
detuning A which represents the energy spacing between 
the ground state | g ) k and the excited state | e ) k of the 



superatom in the rotated frame of reference. The sec- 



ond contribution H, 



coup = ^iJ2k=i a ^' is off-diagonal 
and couples the two levels \e) k and \g) k via the col- 
lective Rabi frequency f2. The remaining part H m t = 

^2k=i SjLfc+i Vkjni^ni^ represents a two-body inter- 
action between Rydberg excitations. More precisely, we 
consider a long-range van der Waals interaction between 
Rydberg excitations, 
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where is the position of the A:-th site in the lattice, i.e., 
the interaction depends on the spatial separation between 
the sites. In this work we assume a repulsive interaction, 
Vk.j > 0, which is common for Rydberg atoms in their 
ns state [25| . In the configuration of a regular ID chain 
with lattice constant a, the distance between two sites k 
and j is simply \j — k\a. Thus, the interaction strength 
reads 
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where we defined V\ := C^/a 6 . 

For the case of regular squares and rectangular lattices 
with spacing a, the N sites are assumed to be arranged 
in M columns and N/M < M rows such that the site 
index k increases as follows: 
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The index k = 1, . . . , N of a site at row = 1, ... , N/M 
and column cj, = 1, . . . , M is given by A: = M(rk — 1) + Ck- 
Conversely, we have = ([(k — 1)/MJ) + 1 and Ck = 
[(fc — 1) mod M] + l, where [^J is the floor function, i.e., 
[x\ is the largest integer not greater than x. The spacing 
between rows and columns is given by the lattice constant 
a and therefore the interaction strength between site k 
and j reads 
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The Hilbert space of Hamiltonian (fTJ) is spanned by the 
tensor product of the single-site Hilbert spaces 
®fcli • The laser contributions -ffdet and H coup are 
local, i.e., they can be conveniently represented within 
%( k ) . The Rydberg interaction, on the other hand, is 
a non-local term and connects different sites. In this 
case the representation is only possible in the full Hilbert 
space T-L. For the investigations in this work it is advan- 
tageous to define the so-called canonical product state 
basis, 
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In this basis the laser detuning part and the interaction 
part of the Hamiltonian are diagonal, whereas the laser 
coupling part of the Hamiltonian is off-diagonal. 

For later purposes it is convenient to introduce the 
excitation number operator, 

N 

^ = E"e fc \ (7) 
fc=l 

which counts the number of Rydberg atoms in the lattice. 
The laser detuning part ifdct then becomes 

H dct = A(N e - N/2). (8) 

In the basis of canonical product states the operator N e is 
diagonal, i.e., the canonical product states | a ) are eigen- 
states of the operator N e with integral eigenvalues N e (a). 
Since we consider in this work only the canonical prod- 
uct states as basis we remark that we usually replace the 
operator N e by its eigenvalue. 

B. Symmetries 

The geometry of the lattice imprints symmetries in the 
Hamiltonian, which we analyze in the following. The 
knowledge about these symmetries can be used to, e.g., 
truncate the state space and to find degenerate states. 

If the Hamiltonian H nlo d commutes with a symme- 
try operator S, the Hamiltonian conserves the symme- 
try, i.e., it couples only states within a certain symmetry 
subspacc. In the matrix representation of the basis of the 
cigcnstatcs of the symmetry operator, this yields a block- 
diagonal form of -ff m od- Correspondingly, when initially 
preparing the system in an eigenstate of the symmetry 
operator, the time evolution of this state is restricted to 
its symmetry subspace. Hence, in such a case it is not 
necessary to consider the full state space. Rather it is suf- 
ficient to restrict the consideration to the corresponding 
subspace. 

Throughout this work we presume the canonical 
ground state | G ) = | ggg . . . g ) as initial state. The ac- 
tual relevant subspace is then determined by the sym- 
metry properties of \G); depending on the number of 
symmetries, this reduces the number of necessary basis 
states and therefore the numerical effort significantly. We 
will give some exemplary numbers in the end of this sec- 
tion, but first we consider the symmetries of different 
lattices and indicate how to find a basis for the necessary 
subspace. 

a. Chain For a linear chain with constant lattice 
spacing and global laser parameters we define a reflec- 
tion operator Y by its action on the Pauli-matrices, 

y ( xWy- 1 =4f +1 " fc ), me{x,y,z}, (9) 

i.e., it exchanges site k with site N + 1 — k. By construc- 
tion we have Y 2 = 1 and the two possible eigenvalues are 
{+!,—!}. Each of the 3 terms in the Hamiltonian i? mo d 
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(a) : 

(b) (c) 




FIG. 1. Symmetry operations for (a) the chain, (b) the square, 
and (c) the rectangular lattice. 



commutes with this operator Y, i.e., [i? mo d,^] = 0. The 
canonical ground state | G ) is obviously an eigenstate 
to the reflection operator with eigenvalue +1, Y\G) = 
+1|G), such that it is sufficient to truncate the Hilbcrt 
space to states which are also cigcnstatcs of Y with eigen- 
value +1, the so-called symmetric subspace. 

When symmetrizing every canonical product state | a ) , 



|a sym ) =Af(l + Y)\a), (10) 



the set of disjoint | a sym ) forms a basis of the sym- 
metric subspace. Af is a normalization constant such 
that ( a sym | a sym ) = 1 holds. The reader may eas- 
ily verify that the symmetrized states are eigenstatcs 
of the reflection operator with eigenvalue +1. The 
symmetrized state | a sym ) has contributions of either 
one or two different canonical product states depend- 
ing on whether the canonical product state is already 
an eigenstate of the symmetry operator Y", e.g., \ct) = 
I egg gee) -> |a sym ) = ^ (\ egg gee) + \eeggge)) or 
| a ) = | egeege ) — > | a sym ) = | egeege ). 

b. Square For quadratic lattices we can define a set 
of symmetry operators P = {1, i?, R 2 , i? 3 , X, D, Y, A} 
which all commute with the Hamiltonian, [i? mo d,'S'] = 
0, VS* G P. These arc: unity (1), (p-fold) clockwise rota- 
tions by 90 degrees (R p ), reflections with respect to the 
horizontal and vertical axis (X, Y), and reflections with 
respect to the diagonals (D, A). The set P represents 
the 8 elements of the symmetry group of the square, the 
dihedral group D4, in the state space. In order to define 
the action of the operators conveniently, we momentarily 
replace the site index k with the tuple (rj,, Ck) according 
to its row 7'^ = 1, . . . , M and column = 1, . . . , M in 



4 



the lattice. The actions of the operators then read 

D_.(r fc ,c fc )D-l _„■(<*, M+l-r k ) 

y ff h, c k ) y-i _ (M+x-ffc, c fc ) 



Aat h ' Ck) A~ l 



(r k , M+l-cO 
^(M+l-Cfe, M+l-r fe ) 



(11) 

(12) 
(13) 
(14) 
(15) 



X" 



where m € {a;, y, z}. By definition we have R A 
Y 2 = D 2 = A 2 = 1. It is important to keep in mind that 
the mentioned operators do not commute in general. In 
Fig. [T]we have indicated the symmetry operators. 

The canonical ground state | G ) is an cigenstatc of all 
8 symmetry operators with eigenvalue +1, 



S\g...g) =+l\g...g), VSeP. 



(16) 



Accordingly, the Hilbcrt space can be truncated to eigen- 
states of all 8 symmetry operators with eigenvalue +1; we 
call this subspacc symmetric as well. 

By means of the group property SP = P, VS e P, one 
finds that S' (EseP S ) = EseP S > V5 " e P - Hence, for 
any canonical product state | a ) the superposition 



(17) 



\S£P 



is an eigenstate to all symmetry operators with eigen- 
value +1, where J\f is again a normalization constant. 
As in the case of the chain, the resulting set of disjoint 
| a sym ) forms a basis of the symmetric subspacc. The 
number of different terms in Eq. (fl"7f is 1, 2, 4 or 8 de- 
pending on whether the canonical product state | a ) is 
already an eigenstate to (multiple) symmetry operators 
in P. These numbers are divisors of the number of group 
elements in P and result from group theoretical consid- 
erations in the context of orbits and stabilizers [26l j . 

c. Rectangle The case of the rectangular lattice can 
be treated analogously to the square lattice. Here, 
we have a set P containing 4 symmetry operators 1, 
Riso, A, Y which represent the symmetry group of the 
rectangle, the dihedral group D 2 - X and Y are horizon- 
tal and vertical reflections, respectively, and i?iso is a 
rotation of 180 degrees. Using again the notation (?*&, c&) 
with r k = 1 . . . N/M and c k = 1 . . . M, we find 



R 



r>, Cfc) -y—l - (N/M+l-r k , c k ) 
y (r k ,c k )y-l _ (r k , iW+l-Cfc) 

(N/M+l-r k , Af+l-Cfc) 



ISO 



n {r k ,c k )y-\ 
°™ 1 180 



(18) 
(19) 
(20) 



where m e {x,y,z}. All four symmetry operators com- 
mute with the Hamiltonian and satisfy R 2 S0 = X 2 = 
Y 2 = 1. The canonical ground state is again an eigen- 
state to all symmetry operators with eigenvalue +1 and 
Eq. (JTZJ) holds with P = {t,R 1S0 ,X,Y}. Accordingly, 
the number of different terms in Eq. (|17p for the rectan- 
gle is 1, 2 or 4. 



TABLE I. Number of basis states of the full Hilbert space 
and of the symmetric subspace of the Hamiltonian. For the 
rectangular lattices we assume a square lattice, but only apply 
the symmetries of a rectangle to truncate the state space. 



Number of sites 


N = 4 


N = 9 


N= 16 


AT = 25 


Full state space 


16 


512 


65536 


~ 34 x 10 6 


Linear chain 


10 


272 


32896 


~ 17 x 10 6 


Square lattice 


6 


102 


8548 


~ 4 x 10 6 


Rectangular lattice 


7 


168 


16576 


~ 8 x 10 6 



The number of basis states of the full Hilbcrt space T-L 
is given by 2 N , N being the number of sites. As outlined 
before, the number of basis states reduces when consid- 
ering the introduced symmetric subspaces. In Tab. U we 
provide specific examples for the number of basis states 
of the symmetric subspaces. For larger lattices, the state 
space is approximately reduced by a factor of 2, 4 and 
8 for chains, rectangles and square lattices, respectively. 
For large square lattices, the majority of symmetrized 
states | a sym ) consist of 8 different summands in Eq. (fl7|) 
since most of the canonical product states | a ) are a pri- 
ori not eigenstates to any symmetry operator. The few 
cases with less terms are negligible and the state space 
reduces by a factor of almost 8. For the chain and the 
rectangle the situation is similar. 



III. SPECTRA 

In this section we focus on the spectra of the symmetric 
subspace of square lattices for a weak laser coupling, i.e., 
|fl| <C Vi,|A|. We also compare the results with the 
case of a linear chain which has been already discussed 
thoroughly in (l9j . In the weak laser coupling regime the 
Hamiltonian is dominated by the diagonal contributions 
-ffdet and -ffint, while the laser coupling term H coup leads 
to a small off-diagonal perturbation. Employing Eq. © , 
the diagonal part of the Hamiltonian reads 



iT diag = A(N e - N/2) + 



N-l N 

E 

i j=k+i 



V kJ n 



(21) 



and its energy eigenvalue for a canonical product state 
| a ) is given by 



E(a) 



Vi 



(N e (a) - N/2) + 



E hlt (a) 
V 1 ' 



(22) 



where E- ini (a) is the eigenvalue of Hi ni , i.e., the inter- 
action energy of the state \a). In Fig. [Ha,c) we show 
the spectrum of Hamiltonian (|2"Tj) , i.e., the A-dependent 
energy eigenvalues, for a linear chain and for a square 
lattice. Obviously, each state appears as a straight line 
in the spectra since E(a) is linear in A. Depending 
on the specific ratio A/Vi we observe points of high 
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degeneracy (crossings) between states. According to 
Eq. (|2"2"j) , two canonical product states | a ) and | (3 ) with 
N e (ct) N e ((3) become degenerate at 



Vi 



N e (a) - N e (J3) 



(23) 



This means that the laser detuning can be used to se- 
lectively compensate for the level shift due to Rydberg 
interactions in order to achieve a degeneracy between a 
pair of states \a) and |/3). If N e (a) = N e {j3) the states 
are either degenerate for all A [for E- m t(a) — Ei nt (f3)] or 
never become degenerate [for Ei nt (a) ^ E- mt (/?)]. In the 
following, we discuss the case of the linear chain and the 
square lattice in more detail. 

a. Chain The case of the linear chain has been dis- 
cussed in [l9[ . In the regime of weak laser coupling it 
is sensible to approximate the interaction energy Ei nt (a) 
by next-neighbor interactions only. The number of neigh- 
boring excitations of a state \a) is denoted by N ee (a), 
e.g., \ a) — \ eeegee) — > N ee (a) = 3. According to Eq. §5§ 
neighboring excitations contribute V\ to the interaction 
energy, thus we have E- mt {a) « N ee {a)V\. Eq. (|2"3"j) then 
reads 



A N ee {a) ~ N ee (P) 
Vi N e (a)-N e (P) ' 



(24) 



Since the numerator and the denominator are integral 
numbers we observe points of degeneracy in the spectra 
at rational A/Vi. The large spacing between degenera- 
cies at integral A/Vi in Fig. [2][a) is a result of the dis- 
cretization of -Ej n t in multiples of V\ . For fixed but inte- 
gral A/Vi all a-dependent terms on the right hand side of 
Eq. (|22p are in the next-neighbor approximation integral 
numbers. As a consequence, also E(a) jV\ is discrete with 
spacing 1. For rational values AjV\ = — |, p,5 6 ^\{0} 
and \p\, \q\ coprime, the spacing in E(a)/Vi in Eq. (|2"2l) is 
reduced to A, which explains the smaller gaps between 
degeneracies visible in Fig. [Ha) for non-integral A/Vi. 

b. Square lattice The situation for the square lat- 
tice is similar to that of the linear chain. However, one 
needs to take into account diagonal neighboring exci- 
tations to approximate the interaction energy. Anal- 
ogously to N ee (a) we introduce the number of diago- 

r /TJi 

nal neighboring excitations Nee («), which have a spa- 
tial separation of \[2a. The interaction energy in next- 
and diagonal- neighbor approximation reads E\-^(a) = 

[N ee (a) + lN e ^ ] (a)]Vi. To follow the arguments of the 
ID case we introduce the effective number of diagonal 
neighboring excitations 

N^%(a):=8N ee (a)+NifKa), (25) 

which is also an integral number. Hence, £"i n t(Q;) = 
\N}^} R (a)Vi 1 i.e., the interaction energy is also dis- 
cretized but on a finer grid with spacing Vi/8. Accord- 
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Laser Shift A/V 1 

FIG. 2. (a) Spectrum of the chain with 16 sites. (b,c) Spec- 
trum of the 4x4 square lattice, where (b) accounts only for 
next- and diagonal-neighbor interactions. For the square lat- 
tice the pattern in the interaction energy occurs on a finer 
grid compared to the linear chain: the spectrum shows large 
energy gaps of |Vi whenever A/Vi is an integral multiple of 
^. Due to long-range interactions that lift the degeneracies of 
energy levels, this pattern is hardly visible for higher energies 
in the square lattice, cf. upper half of subfigure (c). The axes 
of the subfigures have different scales to ensure the visibility 
of the structure. 



ingly Eq. 



reads 



A 



<l^)-A^k/3) 

' Ne(a)-N e ((3) ■ 



(26) 



and we expect again degeneracies between states | a ) and 
I/?) at rational values of A/Vi, cf. Fig. [H[b). Insert- 
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ing A/Vi = — 1|, p,q E ^\{0} and coprime and 



£int(a) = l^llrW in E( l- © wc find 



£(a) 



--(iV e (a) 



Since N e (a) and A^ e e g(a) are integers we expect a spac- 
ing in E(a)/Vi of |t|t. Equation (|27)l describes why we 
observe much smaller energy gaps between degeneracies 
in square lattices and that energy gaps of |Vi appear at 
multiples of 1/8 for A/V u cf. Fig. [2(b). 

1^/21 

In general, N e , N ee , and Nee are restricted by the 
geometry of the lattice and not all combinations are pos- 
sible. As a result we partially observe larger gaps between 
degeneracies in the lower half of Fig. [2(b) than expected 
from our previous argumentation. Specifically, due to the 
small size of the 3x3 lattice, its spectrum shows on the 
one hand the fine structure of the square lattice and on 
the other hand still large gaps at integral A/Vi like in 
the case of the chain, cf. Fig. [3J 

Let us now have a look at the influence of the long- 
range interactions neglected so far. In Fig. [2(b,c) we 
compare the spectra of the 4x4 lattice for next- and 
diagonal neighbor interactions to the one including the 
long-range interactions. Although still visible, the struc- 
ture in the spectrum is less developed, i.e., the gaps di- 
minish, especially in the upper half of subfigure (c). This 
can be explained as follows: The contributions of long- 
range interactions wash out the discretization of the in- 
teraction energy E int (a) = |iV] e ^(ct^Vi. Therefore the 
lines in the spectrum (b) turn partially into a large num- 
ber of narrow lines which are recognizable as thick lines in 
subfigure (c) . Points of high degeneracy turn into regions 
with a high density of states and a large number of degen- 
eracies. As a result also the spacing in E(a) is washed out 
and the gap sizes in the spectrum diminish. As we can 
see in the upper half of subfigure Fig. [5(c), already for the 
4x4 lattice the deviations due to long-range interactions 
are large enough to partly fill even the gaps of |Vi in the 
spectrum. We expect that these effects become stronger 
with increasing lattice size due to an increasing number 
of long-range interactions. The same applies in principle 
also to the linear chain. However, in this case the ef- 
fects of long-range interactions are smaller compared to 
square lattices since the gaps in the chain are larger (up 
to Vi) while the next-to- next- neighbor contributions are 
of order V\/QA. In square lattices we usually have gaps of 
V\ /8 but the next-to-diagonal- neighbor contributions are 
of order Vi/64, as well. In addition, these contributions 
are more numerous in the case of the square lattice. We 
remark that the long-range interactions shift windows in 
the spectrum slightly to more negative values of A/Vi, 
cf. Fig. He). 



N/2) 



(27) 
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FIG. 3. Spectrum of the symmetric subspace of Hamiltonian 
([T]) for the square lattice with 3x3 sites with a laser coupling 
of Q — 0.05Vi . The lower half of the plot shows the projection 
probabilities of the ground state (red, dashed) on a sample of 
(symmetrized) basis states. 



IV. GROUND STATES 

In this section we focus on the laser detuning- 
dependent ground states of the symmetric subspace of 
the system. This is of interest, e.g., for an adiabatic 
sweep through the ground states to selectively prepare 
a desired state [2l| . In Fig. [3] we indicate the energetic 
ground state in the spectrum of the 3x3 lattice with 
a dashed, red line. The associated projection probabili- 
ties Pi, also shown in Fig.[3j have been obtained by exact 
diagonalization and obviously depend on the laser detun- 
ing A/Vi. In general, a numerical diagonalization is not 
feasible due to the exponential growth of the state space. 
We therefore consider the limit of a vanishing laser cou- 
pling, O = 0, which avoids the diagonalization since the 
resulting Hamiltonian 



N 



diag 



N 



fe=l 



JV-1 

E 

k=l j=k+l 



E ^ 



(28) 



is diagonal in the basis of canonical product states. Using 
Hamiltonian (|28[). the projection probabilities in Fig. [3] 
turn into step functions rather than having smooth edges. 
Even in this case the determination of the ground states 
is non-trivial if one refrains from evaluating the complete 
diagonal, which also becomes inefficient with increasing 
lattice size. A more sophisticated method to obtain the 
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ground states of Eq. (|28]l is to solve the problem via 
maximum cuts as we will explain in the next section. 
Prior to that we summarize some general properties of 
the ground states. As in the previous section, the eigen- 
values of -ffdiag are given by Eq. (f22j) . When searching for 
the ground state for fixed A/Vi we can neglect the con- 
stant offset — iVA/(2Vi). Hence, we are interested in the 
canonical product state | a ) which minimizes the energy 



Vi Vi y ' Vi 



(29) 



To achieve this, we divide the canonical product states 
| a ) into groups according to their number of excitations 
N e (a) such that within each group N e is fixed. Thus, for 
each group the term N e {a)A/V\ in Eq. (|29|) is constant 
and minimizing E(a) for fixed laser detuning A implies 
minimizing the interaction energy Ei nt (a). As a result, 
the configuration of minimum interaction energy (MEC) 
for each possible N e = 0, . . . , N is a ground state candi- 
date; vice versa, every ground state is a MEC. Determin- 
ing the true ground state from the N + 1 MECs is then a 
trivial task. Finding the N + 1 MECs, however, is a non- 
trivial optimization problem. We remark that there is in 
general more than one MEC per N e since states are usu- 
ally degenerate for symmetry reasons or might be even 
accidentally degenerate. From Eq. (f2U)) it is evident that 
for A/Vi > the ground state is given by the canonical 
ground state \G) since Ei nt (a) > 0. With decreasing 
A/V\ the number of excitations of the ground state in- 
creases monotonically and for small enough A/Vi <C 
the fully excited lattice \ee. . .e) becomes the ground 
state. See Fig. [3] for an example. 



A. Ground states and max-cut 

In this subsection we show that the evaluation of the 
ground state of Eq. ([28]) for a certain A/Vi is equivalent 
to a combinatorial optimization problem, the maximum 
cut problem (max-cut). A sophisticated code for solving 
max-cut allows then to obtain the ground states for a 
large number of systems. For the formulation as max- 
cut problem it is necessary to transfer our Hamiltonian 
to the Ising-modcl form, 



N N-l N 

Hi s (uj) = ~y] h k Sk - ^2 2J JkjSkSj 

k=l j=k+l 



(30) 



k=l 



Here, the spins take on only two discrete values S k = 
±1 and the spin configuration is indicated by uj. The 
strength of the interaction is given by real weights J k ,j 
and hk is a site-dependent external magnetic field. 

For the transformation to the Ising model form, we 
first rewrite -f/diag in terms of Pauli matrices. In ap- 



pendix [X] we demonstrate that using 



,(fe) 



r (fc) 



t}/2 



Hamiltonian yields 

fc=i \ j^k j 

N-l N 

+ E E 4 V «A k) vi j) + const., (31) 
fc=i j=k+i 

where Vkj = V(\rk — | ) . The desired Ising model form 

(|3"0l can be obtained by translating ai — > Sk while 
omitting the constant offset: 



N-l N 



^=E(t + zE^I^+E E \y^s kSi . 

(32) 



,24 

k=i \ j^k 



k=l j=k+l 



In a physical interpretation, the Rydberg-Rydberg in- 
teraction leads to an additional site-dependent external 
magnetic field in the Ising model. 

If we define a fixed 'ghost spin' Sq = +1 and put Jq j = 
hj , Hamiltonian (|32j) can be written in the compact form 



N-l N 



H ls (u) = - E E J k,3 S kSj, 



(33) 







k=0 j=k+l 


with 






Jo,j = 




1 N \ 






i=i j 




k^J, k,j i v. 



As presented in [27( the ground states of such an Ising 
Hamiltonian can be obtained by means of the concept of 
max-cuts in graph theory. For this, the individual spins 
k in the lattice are represented by vertices (or nodes) 
k in the vertex set V of a graph G. The interaction of 
two spins k,j is represented by the edge (k,j) connecting 
the nodes k and j, its ends, in the graph. The interaction 
strength is taken into account by a real edge weight Ckj — 
— Jkj- The choice Sk € {+1, —1} divides the set of nodes 
V into two (disjoint) subsets V + = {k 6 V\Sk = +1} and 
V- = {fc € V\S k = -1}, respectively. A cut S(V + ) = 
S(V~) is given by the set of edges, where one end is in 
V + and the other in V~ ; the weight of a cut is given by 
the sum of the weights of all its edges. As demonstrated 
in [27| , minimizing Eq. (|33[) is equivalent to finding a cut 



and rearranging the terms of the interaction part of the 



with maximum weight. We add some remarks: 

• The idea given above is valid for a large class of 
systems since the geometry and the interaction 
potential enter only via the interaction potential 
V kd = V(\T k - rj \). 

• Max-cut will only return one configuration u for 
which Hi s (uj) reaches a minimum. However, in gen- 
eral this configuration is degenerate with others and 
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it is not predictable which of them max-cut will re- 
turn. If the degeneracies occur for symmetry rea- 
sons, obtaining the other configurations is possible 
by applying the symmetry operators of the lattice, 
cf. Sec. Ill B[ An additional accidental degeneracy 
(not for symmetry reasons) remains undetected in 
general. 

• The max-cut problem is in general NP-hard [27j 
and therefore also challenging. 

Using max-cut to determine all ground states of the 
system requires to solve the max-cut problem many times 
in a complete parameter range of A/Vi. We scanned 
the parameter range recursively and the recursion was 
stopped by limiting the maximal resolution in A/Vi to 
a certain threshold, e.g., 10 -6 . This means if a MEC 
becomes a ground state for a range in A/Vi which is 
below this threshold, we will in general not detect this 
MEC. The same applies to configurations which never 
become ground state of the system. The latter case ap- 
pears for example in 2D lattices. Nevertheless, for most 
of the below considered systems we can present the com- 
plete set of N + 1 MECs, because we additionally deter- 
mined the MECs with a brute force method. For this the 
complete Hamiltonian diagonal is evaluated which takes 
for 36 sites on a single-core workstation about 2 weeks 
and the computation time roughly doubles with each ad- 
ditional site. In comparison, a corresponding max-cut 
instance can be solved in only 10 minutes, being an im- 
pressive enhancement such that max-cut allowed us to 
determine the MECs of lattices with up to 49 sites. 



B. Ground state sequences 

In the previous section wc introduced our method to 
search for ground states. In this section, we present re- 
sults for selected examples. We start with a brief dis- 
cussion for the linear chain and then turn to the ground 
states of the square lattices. Subsequently we summarize 
some properties wc found in the case of a rectangular 
lattice and show two examples. 

In the following, we will often talk about MECs: as a 
function of the laser detuning A, the ground state is a 
sequence of MECs, each of which is the ground state of 
the system for a given range in A. 



1, Cham 

We start with the discussion of the ground states of 
the linear chain. For the linear chain we are able to cal- 
culate the MECs for up to 36 sites with max-cut and 
verified them by the brute force method (double preci- 
sion: ~ 10 -15 ). The difficulties for the max-cut code to 
find the ground state of larger linear chains stem from 
the very small interaction contribution if the excitations 
are far apart from each other, e.g., the smallest one being 
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iVe = 
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-1.5 -1.0 -0.5 0.0 
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FIG. 4. Illustration of the MECs/ground states of the linear 
chains with TV = 15 (a) and TV = 16 (b). Excitations are 
filled (red/gray), ground state atoms in white. Note, these 
states are usually degenerate for symmetry reasons, (c) and 
(d) show a detail of the spectrum, where the ground state is 
indicated by a red, dashed line. 



Vi/(N — l) 6 . The necessary accuracy to distinguish be- 
tween a ground state and an excited state poses a numer- 
ical challenge and therefore reduces the possible/feasible 
chain sizes. 

In Fig. Q|a,b) we present the ground state MECs for 
the chains with N = 15 and N — 16. The crossover laser 
dctunings where the ground state switches from one MEC 
to another are given in Tabs. UI1 and HTT1 Note that these 
states are usually degenerate with other canonical prod- 
uct states due to the reflection symmetry of the chain; the 
ground states of the symmetric subspace can be obtained 
by employing Eq. (TlU)) . Product states with N e = 1 are 
exceptions because a single excitation can be placed on 
an arbitrary site without changing the interaction energy. 
Thus, the N canonical product states with = 1 are all 
degenerate and the ground state of the symmetric sub- 
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TABLE II. Chain with N = 15: Range for which the MECs 
in Fig. UK a) become a ground state. For example, the MEC 
with N e = 15 becomes a ground state in the range — oo < 
A/Vi < -2.034670. 



A/Vi 


N e 


A/Vi 


N e 


A/Vi 


N e 


— oo 


15 


-1.989622 


9 


-0.000358 


3 


-2.034670 


14 


-1.986998 


8 


-0.000017 


2 


-2.034494 


13 


-0.045188 


7 


0.000000 


1 


-2.033691 


12 


-0.044529 


6 








-2.029601 


11 


-0.017995 


5 


oo 




-2.003177 


10 


-0.002879 


4 






-1.989622 




-0.000358 









N f: = N e = l N e = 2 N e = 3 N e = 4 




iV e = 5 (JV e = 6) (A r „ = 7) jV e = 8 N e = 9 

FIG. 5. Illustration of the MECs/ground states of the square 
with N = 9. Excitations are filled (red/gray), ground state 
atoms in white. Note, these states are usually degenerate for 
symmetry reasons. Parentheses indicate MECs which do not 
appear in the ground state sequence. These have not been 
found by max-cut but by a brute force method. 



TABLE III. Chain with N = 16: Range for which the MECs 
in Fig. |4jb) become a ground state. 



A/Vi 


Ne 


A/Vi 


N e 


A/Vi 


N e 


-co 


16 


-1.989913 


10 


-0.001927 


4 


-2.034673 


15 


-1.988232 


9 


-0.000182 


3 


-2.034555 


14 


-1.016109 


8 


-0.000012 


2 


-2.033875 


13 


-0.045014 


7 


0.000000 


1 


-2.030749 


12 


-0.044340 


6 








-2.016152 


11 


-0.004830 


5 


CO 




-1.989913 




-0.001927 









space is a superposition of symmetrized product states 
with N e = 1. Besides the symmetry- induced degenera- 
cies, in general it might be possible to have states being 
accidentally degenerate with the obtained MECs. 

The MECs in Fig. |U[a,b) seem to trivially separate 
the excitations as far as possible, especially for N e < 
N/2. However, due to the interplay of many long-range 
interactions, it is hard to predict the configuration of 
minimum interaction energy, in particular if the chain 
size increases. Hence, we will not further discuss the 
excitation patterns of the MECs but consider them in 
the following as given. 

Empirically, we find for all investigated linear chains 
that for each number of Rydberg atoms N e there is a 
MEC which becomes the ground state for some laser de- 
tuning A/Vi. This is in general not the case for other lat- 
tices as we will show later. However, most of the ground 
states appear only for very small intervals in A/Vi and 
could be therefore hard to detect in an experimental real- 
ization. In other words, while varying the laser detuning 
we find regions where the ground state changes rapidly 
(clustered ground state crossovers) and regions where the 
ground state does not change in a wide range, cf. dashed 
line in Fig.[4][c,d). For the examples N = 15 and A" = 16, 
most dominant ground states are the MECs with N e = 8 
and N e = 8, 9, respectively and trivially N e = and 
N e = N. This is also visible in Fig. 0Jc,d) and in Tabs. 



HIl and IIIII (bold numbers) . 

In general we observed that for the chain the MECs 
with N e = and N e = \^^-~\ become ground states 
in a wide range in A/Vi, where \x~\ is the ceiling function, 
i.e., \x] is the smallest integer not less than x. To qualita- 
tively understand this, we have a closer look at the energy 
difference between two subsequent MECs | a ) , | (3 ) with 
N e (a) and N e (f3) = N e (a) + 1 excitations, respectively. 
The laser detuning where E(f3) = E(a), i.e., for which 
the ground state changes from the MEC with N e (a) to 
N e ((3), reads according to Eq. (j2"3")l 

AAr e ( a ),jv e ( Q ) + i = - [E int (/3) - E int (a)] , (36) 

where the ground state for A < A7v e (Q),7Ve(Q)+i is repre- 
sented by |/3) and for A > Ajv e ( a ),jv e (a)+i by |a). 

For the case of N = 15 we consider the MECs with 
N e = 7,8,9, 10. The associated interaction energies are 
Eint/Vi rj 0,0,2,4, respectively, approximately given by 
the number of neighboring excitations. Hence, we arrive 
at the crossover detunings A7 i g/Vx ~ 0, Ag^/Vj w —2 
and Agjo/Vi w —2 agreeing with Tab. HIl and the wide 
range in Fig. @|c) for which the MEC with N e = 8 is a 
ground state. For A^ = 16 we find Ei nt /Vi ~ 0, 0, 1, 3 for 
the MECs with N e = 7,8,9,10. This yields A 7 , 8 /Vi « 
0, A 89 /Vi ps -1 and A 9 W /Vi ps -2, cf. Tab. inland 
Fig.HTd). 

In general, if we assume that the excitations of the 
MEC with A^ = A^ G f a chain with an odd number of 
sites N are equidistantly distributed, this excitation pat- 
tern has N ee = 0, i.e., Ei nt /Vi ps 0. Therefore the follow- 
ing MEC has N ee = 2, E- mt /Vi ps 2 and the alternating 
chain is a dominant ground state. Similar arguments are 
possible for the case of an even number of sites. 

2. Square lattice 

In this paragraph we focus on the ground states of 
square lattices. The max-cut code allows us to determine 
the MECs for squares with up to 49 sites; the MECs are 
verified by the brute force method up to 36 sites only due 
to computation time. For the 7x7 lattice, the resolution 
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TABLE IV. Square with N = 9: Range for which the MECs in 
Fig. [S] become a ground state. Bold numbers indicate states 
which are ground state for a wider range in A ('dominant 
ground states') and the lack of a certain N e in the ground 
state sequence is emphasized by underlining the next higher 



A/Vi 


N e 


A/Vi 


N e 


A/Vi 


iVe 


— oo 


9 


-0.5 


4 


-0.001953 


1 


-4.5 


8 


-0.034781 


3 








-2.698417 


5 


-0.029672 


2 


oo 




-0.5 




-0.001953 









TABLE V. Square with N = 16: Range for which the MECs 
in Fig. become a ground state. 



A/Vi 


N e 


A/Vi 


N e 


A/Vi 


N e 


— oo 


16 


-3.042452 


10 


-0.142953 


4 


-4.565203 


15 


-2.019825 


9 


-0.003003 


3 


-4.440203 


14 


-2.019654 


8 


-0.002655 


2 


-3.170405 


13 


-0.533203 


7 


-0.000171 


1 


-3.169405 


12 


-0.411706 


6 








-3.043452 


11 


-0.170292 


5 


oo 




-3.042452 




-0.142953 









IB 



N e = JV. = 1 A 7 . = 2 AL = 3 iV e = 4 AT. = 5 



llllH 



N P = 6 A P = 7 A e = 8 AT. = 9 AL = 10 N r = 11 



12 N e = 13 AT, = 14 N P = 15 AT. 



16 



FIG. 6. Same as in Fig. [5] but for N = 16. 



in A/Vi of max-cut was set to 10 -8 , i.e., if the program 
misses MECs, they become a ground state for a range 
in A/Vi of less than 10~ 8 . In Figs. [5]|9] we give the re- 
sulting MECs and in Tabs. ITvTlVIIII we specify the range 
for which the MECs become the ground state. States 
which are ground state for a range of about IVi in A and 
more ('dominant ground states') are additionally empha- 
sized by bold numbers in the tables. Note that we again 
give only one representative MEC; the ground state in 
the symmetric subspace can be obtained by evaluating 
Eq. (fl~7f (except for N e = 1, where the ground state of the 
symmetric subspace is a superposition of symmetrized 
product states with N e = 1, cf. Sec. IIVB ip . 

The most interesting feature of the investigated 2D lat- 
tices is the possible lack of some N e in the sequence of 
ground states, i.e., not every number of Rydberg atoms 
N e becomes a ground state at some value of A/Vi. In the 
tables, we have emphasized such a lacking by underlin- 
ing the next higher N e ; in the figures they are indicated 
by parentheses. A simple example is the 3x3 lattice. 
Its results are given in Fig. [5] and Tab. IIVI The complete 
ground state sequence from positive to negative detuning 
contains the MECs with N e = 0,1,2,3,4,5,8,9 excita- 
tions. The MECs with JV e = 6, 7 never become a ground 
state. 

To get a qualitative understanding of the omission of 
certain MECs, we consider the energy difference between 



N. = N e = 1 N e = 2 N e = 3 N e = 4 N e = 5 



N e = 6 N R = 7 N e =8 N e =9 N e = 10 (iV e = 11) 



■ 



N e = 12 N e = 13 (N e = 14) (N e = 15) N e = 16 N e = 17 




N e = 18 N e = 19 N e = 20 N e = 21 N e = 22 N e = 23 




N e = 24 N e = 25 



FIG. 7. Same as in Fig. [5] but for AT = 25. 



two MECs | a), |/3) with N e (a) and N e {f3) = N e (a) + n 
excitations (n = 1,2, .. .), respectively. Furthermore we 
assume that | a) is a ground state for some detuning A'. 
If we want to know for which laser detuning A the ground 
state changes from | a ) to | /3 ) we need to fulfill 

> E((3) - E(a) =nA + E int (f3) - E int (a), (37) 



according to Eq. (f2"2"j) . We abbreviate the detuning where 
E(j3) = E(a) by ^-N e (a),N e (i3)- For our example, we 
set now | a } to be the" MEC with N e = 5 of the 3 x 3 
lattice in Fig. which has E int (a)/Vi « 0.6. The 
MEC with N e = 6, denoted momentarily as |/3i), has 
-Eint(/3i)/Vi ~ 3.6, roughly 3Vi more than \ a). Hence, 
| /3i ) does only have a smaller energy eigenvalue than | a ) 
if A/Vi < A 5 , 6 /Vi w -3. The same, A/Vi < A 5 , 7 /Vi w 
-2.9, applies also to the MEC \f3 2 ) with N e = 7 and 
£int(/3 2 )/Vi ~ 6.4. The MEC |/3 3 ) with N e = 8 has in 
total E[ nt ((3 3 )/Vi w 8.7, which means it reduces the en- 
ergy if A/Vi < A 5)8 /Vi w -2.7. Since the MEC |/3 4 ) 
with N e = 9 has i^int(/34)/Vi w 13.2, its crossover with 
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TABLE VI. Square with N = 25: Range for which the MECs 
in Fig. [7] become a ground state. 



A/Vi 


N e 


A/Vi 


N E 


A/Vi 


N e 


— OO 


25 


-3.066433 


17 


-0.046457 


6 


-4.634313 


24 


-3.066189 


ill 


-0.032038 


5 






O RSI SfiQ 
i.UOlOUO 


13 


U.UU 1 OlO 


4 


-4.526993 


22 


-0.570313 


12 


-0.000543 


3 


-4.508212 


21 


-0.504304 


10 


-0.000464 


2 


-4.134313 


20 


-0.504171 


9 


-0.000031 


1 


-3.070340 


19 


-0.071681 


8 








-3.070095 


18 


-0.046625 


7 


oo 




-3.066433 




-0.046457 









TABLE VII. Square 
in Fig. [8] become a £ 



with N = 36: 
round state. 



Range for which the MECs 



A/Vi 


N e 


A/Vi 


N e 


A/Vi 


N e 


— oo 


36 


-3.053621 


23 


-0.145741 


9 


-4.643048 


35 


-3.053438 


22 


-0.033635 


8 


-4.600027 


34 


-3.028838 


21 


-0.030593 


7 


-4.584402 


33 


-3.028813 


20 


-0.018665 


() 


-4.570631 


32 


-2.020436 


19 


-0.006522 


5 


-4.433069 


30 


-2.020428 


18 


-0.003035 


4 


-4.229609 


29 


-0.574484 


17 


-0.000142 


3 


-4.109845 


28 


-0.542782 


15 


-0.000122 


2 


-3.170676 


27 


-0.537173 


li 


-0.000008 


1 


-3.170650 


26 


-0.394663 


12 








-3.070303 


25 


-0.179229 


11 


oo 




-3.070120 


24 


-0.164336 


10 






-3.053621 




-0.145741 
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N e = N e = l N e = 2 N e = 3 A P = 4 



N e = 5 N e = 6 N e = 7 A e = 8 N e = 9 



N e = 10 N e = 11 N e = 12 (N e = 13) N e = 14 



N e = 15 (N e = 16) N e = 17 N e = 18 N e = 19 



N e =20 N e = 21 N e = 22 A P = 23 N e = 24 



A e = 25 N e = 26 A e = 27 A e = 28 N e = 29 



N e = 30 (iV„ = 31) N e = 32 A e = 33 N e = 34 



A e = 35 N e = 36 



| a) does only appear if A/Vi < As^/Vi w —3.2. Hence, 
when decreasing the detuning from A' on, the ground 
state sequence switches from the MEC with N e = 5 di- 
rectly to the MEC with N e = 8 at A/Vi = A 5 , 8 /Vi sa 
—2.7, in agreement with Tab. II VI In other words, with re- 
spect to the MEC with N e = 5, the MEC with N e = 8 has 
the most beneficial ratio between additional interaction 
energy and additional excitations. This emphasizes once 
more that the ground state sequence and the crossover 
detunings crucially depend on the interaction energies of 
the MECs in the considered lattice. 

A similar situation also appears for the other consid- 
ered square lattices with an odd number of sites (N = 25 
and 49, Figs. and EJ for the MEC with N e = \N/2] . 
Qualitatively, this can be explained as follows: In all 
these observed cases the MEC with N e = \N/2] has 
a checkerboard-like structure with excited corners, cf. 
Fig. [71 N e = 13. Adding further excitations to this 
excitation pattern is unfavorable since it increases the 
interaction energy by at least 3 neighboring excitations. 



FIG. 8. Same as in Fig. [5] but for N = 36. 



However, if one rearranges the excitations and forms the 
inverted structure it is possible to add excitations at the 
corners, which increases the interaction energy only by 
2 neighboring excitations. As a result it turns out that 
the MECs with N e = \N/2] + 1, \N/2\ + 2 arc omitted 
and the ground state sequence has a crossover between 
N e = [AT/2] and N e = \N/2] + 3. The lattices show ad- 
ditional omissions of single MECs, but this depends on 
the individual case of the lattice. 

Interestingly, the 4x4 lattice is not missing any MEC. 
This cannot be generalized to all square lattices with an 
even number of sites as the case of the 6x6 lattice demon- 
strates, cf. Fig. [51 
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N e = W e = 1 N e = 2 We = 3 TV,, = 4 JV C = 5 



W e = 13 W e = 14 W e = 15 



N e = 23 W,,. = 24 N e = 25 



W e = 31 W e = 32 N e = 33 



■ 



W e = 38 N e = 39 



iV e = 44 



W e = 45 



N c = 46 



We = 47 



We =48 



W e = 6 We = 8 We = 9 We = 10 W„ = 11 W„ = 12 



W e = 16 We = 20 W„ = 21 




We = 28 We = 29 We = 30 



We = 34 We = 35 We = 36 



N, =40 We = 11 We = 42 



We =49 



FIG. 9. Same as in Fig. [5] but for N = 49. The MECs for 
N e = 7, 17, 18, 19, 22, 26, 27, 43 are not shown since max-cut 
did not yield them as part of the ground state sequence up 
to a resolution in A/Vi of 1CP 8 . Due to the computational 
complexity, we refrained from their determination by brute 
force methods as we did for smaller lattices. 



JV e = N e = 1 N e = 2 W e = 3 N e = 4 



m 



W e = 5 N e = 6 N e = 7 N e = i 



N e = 9 



I I I I I I I I I I I I I I I I 



N,, = 10 We = 11 (W e = 12) (A r e = 13) N e = 14 

MS EH 



N e = 15 N e = 16 W P = 17 W e = 18 W e = 19 



We =20 We = 21 



FIG. 10. Same as in Fig. [5] but for a rectangle N = 3 x 7. 



TABLE IX. Rectangle with N = 21 = 3 x 7: Range for which 
the MECs in Fig. 1101 become a ground state. 



A/1/! 



— oo 
-4.569993 
-4.548927 
-4.491230 
-3.085181 
-3.083228 
-3.053931 
-3.051978 



21 

20 
19 
18 
17 
16 
15 



N e 



-3.051978 
-2.689414 
-0.535250 
-0.502358 
-0.502114 
-0.067022 
-0.048241 
-0.031259 



14 
11 

10 
9 
8 
7 
6 



AM 



Ne 



-0.031259 
-0.016634 
-0.005432 
-0.000916 
-0.000016 


oo 



3. Rectangular lattices 



TABLE VIII 
MECs in Fig. 



Square with N = 
[9] become a ground 



49: Range for which the 
state. 



AM 



-co 
-4.652127 
-4.635938 
-4.618859 
-4.606706 
-4.563284 
-4.553749 
-4.506435 
-4.506287 
-4.140838 
-4.138884 
-4.135408 
-4.081814 
-3.070846 
-3.070831 



N 



49 

48 
47 
46 
45 
44 
42 
41 
40 
39 
38 
37 
36 
35 



AM 



N c 



-3.070831 
-3.070685 
-3.070492 
-3.053305 
-3.053112 
-3.052966 
-3.052950 
-2.681913 
-0.578998 
-0.570054 
-0.546426 
-0.541998 
-0.493574 
-0.089938 
-0.071157 



34 
33 
32 
31 
30 
29 
28 
25 
24 
23 
21 
20 
16 
15 



AM 



-0.071157 
-0.059398 
-0.057176 
-0.047315 
-0.041627 
-0.033062 
-0.006277 
-0.004117 
-0.002765 
-0.000686 
-0.000048 
-0.000041 
-0.000003 



oc 



_N_ 

14 

13 

12 

11 

10 

9 

8 

6 

5 

4 

3 

2 

1 





We end the discussion of ground state sequences with 
a brief outlook on rectangular lattices. Depending on 
the geometry, we were able to evaluate the ground states 
of rectangular lattices with up to 48 sites by means of 
max-cut. 

Due to the limited scope of this work we restrict our- 
selves to present only two examples: iV = 21 = 3x7 and 
N = 24 = 3 x 8. In the first case we have an odd num- 
ber of sites. The results are summarized in Fig. [TU] and 
Tab. lIXI Similar to square lattices with an odd number of 
sites we find that the MECs with N e = 12 and N e = 13 
are omitted in the ground state sequence. The explana- 
tion is analogous to the discussion of the 3x3 lattice: The 
MEC with N e = 11 forms a checkerboard-like pattern 
while the excitations in the MEC with = 14 resembles 
the inverted pattern with excited corners. The associated 
interaction energies lead to a omission of the intermediate 
MECs. We remark that also the other investigated rect- 
angular lattice with 3 rows and an odd number of sites 
systematically skip the MECs with N e = + 1 and 
A^ e = [-y] +2, since they undergo similar arrangements 
of the excitations as the 3x7 lattice. 
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N c — N e = 1 iVc = 2 N e = 3 JV e = 4 



7V e = 6 JV e = 7 W e = 8 N e = 9 



JV. = 10 iV e = 11 iV e = 12 N, = 13 JV e = 14 



^ffl ^M MM MM 



(N e = 15) jV e = 16 (K. = 17) jV e = 18 N e = 19 



JV e = 20 


Ne = 


21 N e = 22 


N e = 23 N e 


= 24 


FIG. 11. 


Same as in Fig. [5] but for a rectangle N = 


3x8. 


TABLE X 


Rectangle with N = 


24 = 3 x 8: Range for which 


the MECs 


in Fig. 


TT] become a g 


round state. 




A/Vi 


N e 


A/Vi 


N e 


A/Vi 


N e 


— oo 


24 


-2.901020 


14 


-0.046089 


6 


-4.570645 


23 


-2.268020 


13 


-0.017949 


5 


-4.562885 


22 


-1.783945 


12 


-0.011288 


4 


-4.508958 


21 


-0.560548 


11 


-0.003907 


3 


-3.536517 


20 


-0.521419 


10 


-0.000582 


2 


-3.083073 


19 


-0.491402 


9 


-0.000007 


1 


-3.078055 


18 


-0.277148 


8 








-3.058263 


16 


-0.046921 


7 


oo 




-2.901020 




-0.046089 









In the case of 24 sites we observe that the MECs with 
N e = 15 and N e = 17 are omitted, cf. Fig.HUand Tab.|X] 
It turns out that the MEC N e = + 3] is also skipped 
in the other investigated rectangular lattices with 3 rows 
and an even number of sites. The reason for this are ar- 
rangements of the excitations comparable to those of the 
lattice with 24 sites. In summary, seemingly the excita- 
tion patterns of the MECs of all rectangular lattices with 
3 rows follow comparable arrangements of the excitations 
with increasing N e . This might be due to the small size 
of the lattice. 



V. SUMMARY 

We explored the symmetries and spectral properties 
of two-dimensional, laser-driven lattices with a main fo- 
cus on square lattices. Exploiting the geometry induced 
symmetries allowed us to truncate the state space when 
initializing the system in the canonical ground state 
| G ) = \gg . . . g), which facilitates the numerical treat- 
ment. 

An investigation of the spectra in the limit of van- 
ishing laser coupling of square lattices showed that the 



spectra of small square lattices show similar features as 
the linear chain. Specifically, we observe in the spectra 
also points / regions of high degeneracy and large energetic 
gaps between them. However these features appear on a 
scale reduced by a factor of 8 which is due to the diagonal- 
neighbor contributions. Because of the corrections orig- 
inating from long-range contributions, these properties 
are less pronounced with increasing lattice size. 

For the investigation of the laser detuning-dependent 
ground states, we have demonstrated that in the limit of 
vanishing laser coupling the model Hamiltonian can be 
written in the form of an Ising-Hamiltonian with a site- 
dependent external magnetic field. The ground states of 
the obtained Ising-Hamiltonian have been evaluated for 
various one- and two-dimensional lattices by using ele- 
ments of graph theory. The method used in this work, 
however, is in principle applicable to arbitrary lattice ge- 
ometries. In contrast to the case of the linear chain, 
where there is a state for every number of excitations 
N e which becomes a ground state at some laser detun- 
ing, we observed that this is in general not the case for 
two-dimensional lattices. 

We emphasize that the results presented here cannot 
be implied to occur in a similar way for larger lattices or 
other interaction potentials since the optimization of the 
interaction energy crucially depends on the lattice and 
interaction potential. 
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Appendix A: Rewriting the Hamiltonian diagonal 



In this appendix we demonstrate that Eq. 



can be 



rewritten by only using a z -matrices as in Eq. pip . 

For the interaction term H lnt 
EpEf=, +1 V k ^n^ in Eq. g 
ni = J[<Ta +1]. Neglecting momentarily the occurring 

constant energy offset C = \ YJk=i Yl!j=k+i V K. 
interaction term _ffi nt reads 



we use 



the 



N-l N 



4 ^ 

k=l j=k+l 



1 N-l N 

k=l j=k+l 



VkM k 



(Al) 



The first term has already the aspired form and we 
rearrange the individual parts of the sum of the second 
part with respect to ai^ . When defining a number r& 
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such that 

JV-l JV 



JV 



E E M^+^HE^ (A2) 
fe=i fe=i 

one finds 

JV fc-1 JV 

r k = e l - • E l > E u - ; -- 1 • • • v - 



j=fc+i j=i i=i 



(A3) 



Please note, the last equality holds only if the interac- 
tion potential satisfies Vk,j = Vj.k, as is the case for our 
system. We demonstrate briefly the equality in Eq. (|A2[) 
by inserting Eq. (|A3|) : 

JV JV / JV fc-1 \ 

E ^ ] = E E + E *s.» -i fc) ( A4 ) 

fc=l k=l \ K j = k+l j = l J 

JV-l JV JV fe-1 

= E E ^^+EEW fc) . 



fc=l j=fc+l 



fc=2 j = l 



(A5) 



Note, we dropped those indices in the last line where 
empty sums occur. The first term is as desired and we 
remain with showing 



JV-l JV 



E E ^=EE^. 

k=l j=k+l k=2 j=l 



(A6) 



which can be done by employing induction. Starting with 
N = 2, Eq. (|A"6j) becomes 



V h2 a^ = V h2 a^. (A7) 
Under the presumption that Eq. (| A6|) holds for jV sites 



we can show that it also holds for N + 1 sites: 

(JV+1)-1 JV+l 

E E vvoP = 
fc=i j=fc+i 

JV / JV 

= E E %^co+% lJV+10 -^+ i ) 

fe=i y j=fe+i 

JV-l JV JV 

= E E v k ,A 3) +Y, v ^-i N+1) 

k=l j=k+l k=l 

fe=2 J=l j=l 
JV+l fe-1 

= EEvi fc) - (as) 

k=2 j=l 

Thus, recalling the constant energy offset C, in summary 
one can write -ff; n t as follows 

^=Ei(e^1 ^ 

k=l \&k J 

JV-l JV 

+ E E VuoMoM+C. (A9) 



fc=i j=k+i 



Together with the laser detuning part -ffdet = 
Y X^fcLi i the diagonal Hamiltonian in Eq. (p8|) is 
equivalent to 



JV-l 



. 2 4 

JV-l JV 



lE E (A10) 



/c=i j=k+i 
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